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ABSTRACT 

We present Suzaku observations of three iron low- ionization broad absorption 
line quasars (FeLoBALs). We detect J1723+5553 (3a) in the observed 2-10 keV 
band, and constrain its intrinsic Ah column density to Ah > 6 x 10 23 cm -2 by 
modeling its X-ray hardness ratio. We study the broadband spectral index, a ox , 
between the X-ray and UV bands by combining the X-ray measurements and the 
UV flux extrapolated from 2MASS magnitudes, assuming a range of intrinsic 
column densities, and then comparing the a ox values for the three FeLoBALs 
with those from a large sample of normal quasars. We find that the FeLoBALs 
are consistent with the spectral energy distribution (SED) of normal quasars 
if the intrinsic iVn column densities are Ah > 7 x 10 23 cm" 2 for J0943+5417, 
A" H > 2 x 10 24 cm" 2 for J1352+4293, and 6 x 10 23 < A" H < 3 x 10 24 cm" 2 for 
J1723+5553. At these large intrinsic column densities, the optical depth from 
Thomson scattering can reach ~ 6, which will significantly modulate the UV 
flux. Our results suggest that the X-ray absorbing material could be located at 
a different place from the UV absorbing wind, likely between the X-ray and UV 
emitting regions. We find a significant kinetic feedback efficiency for FeLoBALs, 
indicating that the outflows are an important feedback mechanism in quasars. 

Subject headings: (galaxies:) quasars: absorption lines — (galaxies:) quasars: 
general — (galaxies:) quasars: individual (SDSS J094317.59+541705.1, SDSS 
J135246.37+423923.5, SDSS J172341. 08+555340.5) — X-rays: galaxies 



1 Homer L. Dodge Department of Physics & Astronomy, University of Oklahoma, Norman, OK 73019, 
USA, morabito@nhn.ou.edu 

2 Department of Astronomy, University of Virginia, Charlottesville, VA 22904, USA 

3 Max-Planck-Institut fiir Astrophysik, Karl-Schwarzschild-Str. 1, D-85748, Garching, Germany 



- 2 - 



Introduction 



Quasars and active galactic nuclei (AGN) are ubiquitous phenomena wherever we look 
in the universe. The subset of quasars known as broad absorption line quasars (BALQSOs) 
are characterized by blue-shifted, rest-fra me UV absorption trou ghs due to gas outflow. 



Even before the first definitive survey in IWeymann et al.l (ll99ll ). BALQSOs provided a 



unique glimpse into the central engines of AGN, where the nature of the intrinsic absorption 
has important implications for the structu re of the central engin e. To classify a quasar 
as a BALQSO, the traditional definition of IWeymann et al.l (119911 ) requires t he absorption 



trough s to be at least 2,000 km s 1 wide. Recent studies have relaxed this, e.g.. lTrump et al. 
( 120061 ). to include quasars with absorption troughs over 1,000 km s _1 in width. Regardless 



of which definition is used, BALQSOs are further divided into those that exhibit absorption 
troughs from only high-ionization species (HiBALs) and those with low-ionization species 
(LoBALs). Most BALQSOs are H i BALs, and most LoBALs also have high-ionization troughs 
in their spectra ( IWeymann et al.lll991t iTrump et al.ll2006l ). LoBALs are further classified 
by the low-ionization species they exhibit: in this paper, the objects surveyed have strong 



iron absorption, and are therefore iron LoBALs (FeLoBALs). T his type o 



LoBALs and 


3ALQSOs make up -4-7% ( 


Dai et al. 


2010b) and -20-40% ( 


Dai et al. 


2008 


Shankar et al. 


2008a: 


Gangulv & Brotherton 


2008: 


^nigee et al. 


2008; 


Maddox et al. 


2008 



BALQSO is 



2010b). while 



Allen et al.lboioh . respectively, depending on the threshold of the trough width. Focusing on 



FeLoBALs is meaningful because of their unique characteristics and the bearing they have 
on the broader picture: either as an evolutionary link in the normal AGN lifetime, or as a 
geometric interpretation of a uniform class of objects. 

Besides the low-ionization absorption lines, the optical continua for LoBALs are more 
reddened than HiBALs and normal QSOs, suggestive of stronge r dust absorption, perhap s 
from large quantities of dust in the vicinity of the central engine (jSprayberrv fc Foltzlll9 92). 



LoBALs in particular are mor e apparently X-ray weak than normal BALQSOs (I Green et al. 
200 ll : iGallagher et al.l l2002al ). There is a large portion, ~80%, not detected in X-rays 



along with lower than expected values for the UV to X-ray lumino sity ratio, a™ ]] 
lower than expe c ted by ~ —0.9), that leads most previous studies, e.g. iGreen et al.l ( 120011 ): 



Gallagher et al.l (2002al). to conclude that there is extreme X-ray absorption in LoBALs. 
Streblyanska et al.l ( 120101 ). however, found that LoBALs in general have lower column den- 
sities (Nh < 10 22 cm -2 ) than HiBALs. This conclusion may have been influenced by their 
X-ray bright sample, required for sufficient quality to perform X-ray spectral analyses. 



1 The broadband spectral index, 2500A to 2 keV point-to-point power-law slope. 



-3 - 



Low- ionization BALQSOs with broad iron absorption lines (FeLoBALs) represent an 
extreme category. They have the reddest continuum spectra, complex UV spectra, and 
possibly largest X-ray absorption column densities. There are a handful of FeLoBALs al- 
ready studied in X-rays, with only a few detections. Almost all were observed using the 
sensitive X-ray telescope, Chandra. The ne arby FeLoBAL Mrk 231 (z = 0.042) was a lso 
detected by XMM-N ewton and BeppoSAX ([Turner fc Kraemerll2003l ; iBraito et al.l 120041 ) in 
addition to Chandra ([Gallagher et al.ll2002bl ). which allowed for the most in-depth studies of 
an FeLoBAL to date. All three studies find the 0.4-10 keV X-ray emission to be mostly re- 
flected and scattered into our line of sight, and the direct X-ray emi ssion is mostly absorbed. 
Using the PDS detector on board BeppoSAX, IBraito et al.l ( 120041 ) detect the direct X-ray 
emission of Mrk 231 at 3cr in the 15-60 keV band, which allows the authors to constrain 
an intrinsic column density Nm ~ 2 x 10 24 cm -2 . Another recent study focused on two 
FeLoBALs is iRogerson et al.l (120 111 ), where the authors calculat ed a m upper limits for both 



objects, comparing them with a large sample of normal AGN (jSteffen et al.l 120061 ) . Other 
studies focused on subsets of quasars, but had FeLoBALs containe d within their scope . 
Quite a few of these focused on radio-detected objects. In one study, lUrrutia et al.l (120051 ) 
selected 12 luminou s red quasars, with the criteria that they were also detected by FIRST 
( iBecker et al.lll9951 ). and used standard aperture photometry (for details see § 12.2ft to conduct 
X-ray analyses. They found that all of their objects are X-ray absorbed to some degree; in 
addition, they found a slightly steeper spectral slope than that for normal quasars, with an 
unweighted mean of T = 2.2 ± 0.4 and a high-energy region weighted mean of T — 2.1 ± 0.5. 
Brotherton et al.l ( 120051 ) selected five radio-loud BALQSOs, three of which are LoBALs, and 
two of those LoBALs exhibited i ron features. The t hree LoBALs were found to exhibit the 
steepest a ox values in the study. iMiller et al.l (120091 ) also selected 21 radio-loud BALQSOs, 
one of which is classified as an FeLoBAL. The authors proposed that the X-ray properties 
are explained by a model that consists of X-ray emission in the disk/corona and linked to the 
radio-jet. Almost all of these studies calculated intrinsic column density for these objects, 
most of which are in the 10 22 - 10 24 cm -2 range. These will be discussed in detail in § 15.21 In 
this study, we will use our X-ray detections and non-detections to constrain X-ray absorption 
column densities. One of our objects is the first X-ray detection (at the 3a level) by Suzaku 
of an FeLoBAL. 

It is not understood why FeLoBALs possess unique physical characteristi cs; evolution 
and geometry are the contending theories. A geometric interpretation (e.g., lElvis 1 120001 ) 
proposes that the broad absorption line region of a BALQSO simply arises from looking 
straight down a narrow outflow, while normal AGN and quasars are the same object viewed 
from different angles. A recent discovery of a LoBAL transitioning from an FeLoB AL within 



a few years rest frame favors special lines of sight, because of the short time scale (IHall et al. 
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20101 ). On the other hand, iLipari et al.l (120091 ) proposed that BALQSOs, and FeLoBALs 
especially, are early progenitors of normal AGN and quasars, at the very end of an extreme 
starburst period and in the process of blowing out material from type II SNe that is Fe II 
rich. This evolutionary t heory is partially supported by the higher fraction of L oBALs at 
high IR luminosities (e.g., Farrah et al. 2007 ; Urrutia et aDl2009 ; Dai et al. 2010b ). A recent 
Spitzer spectral survey of six FeLoBALs ( Farrah et al. 201ol ) showed significant signatures 
of dust and PAH emissions. Spectral modelin gs in the UV regim e also support a large 
covering fraction of the wind for FeLoBALs (e.g. JCasebeer et al.ll2008l ). Other evidence, such 
as the decreasi ng fraction as a function of radio lum inosity, is consistent with a geometric 
interpretation ( IShankar et al.ll2008at iDai et al.ll2010bf ). Therefore, the population of LoBALs 



and F eLoBALs could also be a combination of both evolution and geometry ( IDai et al. 
2010bh . 



Another interesting question is how FeLoBALs contribute to the quasar feedback pro- 
cess. AGN feedback is widely used in galaxy formation models (e.g., Granato et al. 2004; 
Scannapieco & Oh 2004; Hopkins et al. 2005; Shankar et al. 2008b) to explain the co-evolution 
of AGN an d their host galaxies. Feedback can reprod uce such phenomena a s the M-a re- 
lation (e.g., iGraham et al.ll201ll ; iGiiltekin et al.l 120091 ). star formation rates (ISilk fe Nusser 
2010 ). and even the shortfall in the halo baryon fraction ( Silk fe Nusser 201dh . Whether 
FeLoBALs are an evolutionary stage or a geometric interpretation of AGN, studying their 
feedback efficiency allows us to add to our limited picture of their physical characteristics. 

In this paper we present new Suzaku observations of three FeLoBALs, SDSS objects 
J094317.59+541705.1, J135246.37+423923.5, and J172341. 08+555340. 5 (hereafter J0943+5417, 
J1352+4239, and J 1723+5553, respectively). The targets are classified as FeLoBALs by 
Trump et al.l (120061 ). and are t he brightest (K q < 14.4 mag) of the 2 Micron All Sky Survey 
(2MASS) selected BALQSOs (IDai et al.l \2Q0& \ . J0943+5417 is a radio-intermediate quasar 
detected with a peak flux of 1.71 mJy/beam and RMS noise of 0.145 mJy/beam by Faint Im- 
ages of the Radio Sky at Twenty- Centimeters (FIR ST). None of these o bjects are detected 
at the flux limit in the NRAO/VLA Sky Survey JCondon et alJ Il998h . J0943+5317 and 
J 1352+4239 are typical FeLoBALs with extremely complex rest frame UV spectra, whereas 
J0943+5417 has absorption troughs so wide that they overlap and completely suppress the 
UV continuum. We extrapolate our UV data from 2MASS, which samples the rest frame 
optical bands, and then correct to the rest frame UV. The complexity of the observed op- 
tical/rest frame UV spectra makes it extremely difficult to model the continuum. This is 
addressed in § 3. Since the targets all have redshifts z ~ 2, the XIS on-board Suzaku can 
probe the rest frame energy up to 30 keV, where the rest frame 6-30 keV (observed 2-10 
keV) emission is only moderately affected by absorption even if Ah ~ 10 24 cm~ 2 , allowing 
us to detect the FeLoBALs if they are Compton thin. We study the relationship between 
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X-ray and UV luminosities for these three FeLoBALs, and compare our values to a large 
sample of optically selected normal quasars. Our Suzaku observations and data reduction 
are described in §21 followed by calculations of 2500A luminosities in §3. Section 4 presents 
results, followed by discussion in §5. Throughout the paper we assume a cosmology of H =70 
km s" 1 Mpc" 1 , Q M =0.3, Q A =0.7, and k=0. 



Suzaku Observations and Data Reduction 
2.1. Suzaku Observations 



We observed the three FeLoBALs in 2009 with the Suzaku X-ray observatory ( IMitsuda et al 



20071 ). with exposure times ranging between 32-36 ks for each object. Table [TJ lists the ob- 
servation log. During each observation, the standard data mode was used with the pointing 
observation mode. The three working X-ray Imaging Spectrometer (XIS) cameras (XISO, 
XIS1, XIS3, Koyama et al. 2007) were operated in both the 3x3 and 5x5 photon count- 
ing modes with a minimum time resolution of 8 seconds. Suzaku operates in a low-earth 
orbit, 550 km above the earth, completing one period in 96 minutes. The background in 
our observations is low and stable, without flares, consistent with most Suzaku observations 



( lYamaguchi et al.ll2006l ). 



2.2. Suzaku Data Reduction 

Starting with the cleaned event files, we visually inspect the data using DS^ to check for 
any anomalies and ensure that our data extraction is completely contained within the Suzaku 



2 An astronomical imaging and data visualization application ( Jove fc Mandell 2003 ) 



Table 1. The Suzaku Observation Log 





Observation 


Observation 


3x3 mode 


5x5 mode 


Object 


Start (UT) 


Stop (UT) 


exposure (s) 


exposure (s) 


J0943+5417 


14:08:03, May 24, 2009 


08:20:12, May 25 


28,210.81 


6,015.56 


J1352+4239 


23:54:09, Jun 2, 2009 


18:00:19, Jun 3 


23,792.67 


8,289.84 


J1723+5553 


08:49:36, Jun 4, 2009 


19:33:19, Jun 4 


35,944.81 





Note. - 



Only data for the 3x3 mode exists for J1723+5553. 
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field. Using the standard aperture photometry, the half-power diameter region around each 
source is obtained using the Suzaku Technical DescnptioM, 

along with a corresponding 

background annulus. The background annuli are chosen individually for each mode (3x3, 
5x5) to maximize the background area while keeping the inner radius large enough not to be 
contaminated by the source and the outer radius still contained within the field. The effective 
area for the background annulus is measured from the mean of the inner and outer radii of 
the background region. We then use the High Energy Astrophysics Science Archive Research 
Center (HEASARC) tool XSELECT in conjunction with these regions to filter the energy 
to the observed full (0.2-10 keV), soft (0.2-2 keV), and hard (2-10 keV) bands, extracting 
the counts and exposure times. We first combine the 3x3 and 5x5 mode data, then measure 
the source and background counts for the three CCDs, XISO, XIS1, and XIS3, separately. 
We note that XIS1 is a back-illuminated instrument, which has a higher background and 
sensitivity at lower energy ranges. The count rates and uncertainties for XISO and XIS3, 
which are both front illuminated (FI) instruments, are found separately and then combined 
using the Bayesian estimation method to give a total FI count rate and error for the count 
rate. The final count rates and uncertainties are listed in Table |2j 

We first examine the full band, and do not detect the sources; then we examine the 
soft band, with the same result. In the hard band, we do not detect the sources in most 
situations. We marginally detect SDSS J1723+5553 (2.6a and 1.6a) in both of the FI CCDs 
in the observed 2-10 keV band. Combining the data from the two FI CCDs, we detect 
J1723+5553 at the 3cr significant level (Table The FI CCDs have comparable effective 
areas, and lower background rates at 2-10 keV than the BI CCD, which means they are 
more efficient at detecting faint sources in this energy band. Thus the non-detection of the 
source in the BI CCD (XIS1) is expected, given the low X-ray flux. Therefore, we conclude 
that the detection of J1723+5553 in the 2-10 keV band is real. We do not detect the other 
objects in any of the detectors. 

Finally, we use the HEASARC tool PIMMS to convert the count rates into fluxes. Since 
the intrinsic column density of the quasars is unknown, we use a range of 10 20 < Nh < 10 25 
cm" 2 to calculate the un absorbed, intrinsic fluxes. For the galactic column density, the 



Dickey k, Lockmanl (119901 ) average N H values are used (1.3 x 10 20 cm 2 , 1.21 x 10 20 cm 2 , 
and 3.08 x 10 20 cm" 2 , for J09843+5417, J1352+4239, and J1723+5553, respectively). The 
photon index is taken to be T=1.9 (e.g., Reeves & Turner 2000; Dai et al. 2004; Saez et al. 
2008), equivalent to an energy index a=0.9. The unabsorbed fluxes are calculated from the 
BI and FI rates separately using PIMMS and then the upper limits are combined using the 



3 Available at 



http://heasarc.gsfc.nasa.gov/docs/suzaku/prop_tools/suzaku_td/suzaku_td.html 
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Bayesian estimation method. We use just the FI rates in calculations for J1723+5553. Next, 
we estimate the monochromatic 2 keV (rest frame) luminosity using the unabsorbed fluxes 
(or flux limits) in the observed 0.2-2 and 2-10 keV bands. Using the power-law model, the 
flux, F, is given by, 

F= / hv- a dv, (1) 

where /q is a constant found analytically for the 2-10 and 0.2-2 keV ranges. 

From the simple relation in Equation (U we interpolate the value for the monochromatic 
flux density at the observed energy Ex = 2 keV/(l + z) using, 

U x = fa?- ( 2 ) 
We find the 2 keV rest frame luminosity (hkev) from the corresponding observed fluxes using, 

lv,rest = (Y^^ ^^Li ( 3 ) 

where u, obs = u/(l + z) and D L is the luminosity distance. Since f Uj0 bg is calculated for 
a range of intrinsic absorption, we obtain a range of monochromatic luminosities at 2 keV 
(Table E]). 



3. The 2500A Luminosity 



We calculate the rest frame 2500 A luminosity, Z250O) by extrapolating the observed 
2MASS magnit udes (probing r est fra me 4000-7000A for z ~ 2 objects) using the mean 
quasar SED of iRichards et al.l (120061 ). Since all three objects are at z ~ 2, the 2MASS 
bands sample the rest-frame optical bands of the FeLoBALs, where there is less spectral 



Table 2. Count Rates and Uncertainties for Each Object 





Observed 0.2-2 keV 


Observed 2- 


-10 keV 


Object 


FI count rate 


BI count rate 


FI count rate 


BI count rate 


J0943+5417 
J1352+4239 
J1723+5553 


< 3.0 x 10~ 4 

< 2.1 x 10~ 4 

< 2.3 x 10~ 4 


< 6.7 x 10~ 4 

< 9.0 x 10" 4 

< 6.9 x 10~ 4 


< 3.3 X 10~ 4 

< 3.9 x 10" 4 
4.4 ± 1.4 X 10~ 4 


< 8.5 x 10~ 4 

< 9.9 x 10~ 4 

< 9.0 x 10~ 4 



Note. — The 3cr upper limits are reported here, and the 3cr hardband detection of 
J1723+5553. The count rates are within the half-power diameter region of Suzaku, and are 
scaled to full power after the conversion to fluxes. The units are cnt s~ 1 . 
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complexity than the rest-frame UV bands. The Sloan Digital Sky Survey (SDSS) magni- 
tudes sample the rest frame UV of the FeLoBALs, where it is extremely difficult to model 
the continuum due to the presence of complicated absorption troughs and dust extinction. 
Although the observed 2MASS wavelengths are further from the rest frame UV than the 
observed SDSS wavelengths, we do not have to rely on a difficult to estimate model of the 
continuum. The observed optical spectra of J0943+5417 is almost completely suppressed 
(see Figured]), which demonstrates why we use the 2MASS data that circumvents reliance 
on complex corrections. 

We first correct for Galactic extinction caused by the Milky Way for the 2MASS mag- 
nitudes. This extinction is given for the tt-band in the SDSS quasar catalog for each object 
(Schneider et al. 2007). From that value, we obtain the extinction in the 2MASS bands 
using the Milky Way extinction curve, 

^ ( E{\-V) \ 1 

As ~ \e{b=v) + Rv ) rnv (4) 



where the values of color excesses, E(X — V), and Ry = 3.08 are taken from iPeil ( 119921 ). 
After correcting for Galactic extinction, we calculate the observed flux densities at the cen- 
tral J, H, and K s wavelengths, and obtain the rest frame monochromatic luminosities at 
optical bands (at ~6,900A) using Equation [3j We correct for the intrinsic dust extinction 
in these quasar s assu ming E(B — V) — 0.077 with an SMC-like extinction curve, found by 



Reichard et al.l ( 120031 ) to be a good approximation for LoBALs. Since FeLoBALs have higher 
dust extinctions than LoBALs, we possibly under-correct the intrinsic dust extinction for our 
targets. Finall y, we extrapolate the rest-frame optical luminosities to Z 2 5oo using the mean 



quasar SED of iRichards et al.l (120061 ). The parameters and intermediate results used in all 
these calculations are listed in Table HI The results obtained from J, H, and K s bands are 
consistent within 20% for J0943+5417 and J1352+4239, suggesting no significant complexity 
in the corresponding spectra. For J1723+5553, the 2500A luminosities extrapolated from J 
and H bands are consistent within 20%; however, the value extrapolated from the K s band 
is twice as large. It is possible that the observed K s band is contaminated by a strong H a 
emission line. For all three objects, we use the average of the Z 2 5oo values extrapolated from 
the J, H, and K bands, in our following calculations. 



4. Results 

After calculating the intrinsic monochromatic luminosities at 2500A and 2 keV rest 
frame, we measure the broadband spectral index between the UV and X-ray bands, a ox = 
0.3838 x log (/ 2 fcev/^25oo)) an d then compare our a ox values to those of the large sample of 



Table 3. 2 keV rest frame Monochromatic flux and luminosity. 



J0943+5417 J1352+4239 J1723+5553 





Intrinsic 




















N H 


fx 




lx 


fx 




lx 


fx 


lx 




10 20 


< 1.4 x 10" 


32 


< 4.8 x 10 26 


< 1.0 x 10- 


ss 


< 2.7 x 10 26 


< 9.3 x 10 -33 


< 2.7 x 10 26 


Extrapolated 


10 21 


< 1.5 x 10" 


32 


< 4.9 x 10 26 


< 1.0 x 10- 


-32 


< 2.7 x 10 26 


< 9.4 x 10 -33 


< 2.8 x 10 26 


from 


10 22 


< 1.6 x 10~ 


32 


< 5.3 x 10 26 


< 1.1 x 10- 


-32 


< 3.0 x 10 26 


< 1.0 x 10~ 32 


< 3.1 x 10 26 


0.2-2 keV 


10 23 


< 2.9 x 10- 


32 


< 9.7 x 10 26 


< 2.2 x 10- 


-32 


< 6.0 x 10 26 


< 2.0 x 10~ 32 


< 6.0 x 10 26 


Observed 


10 24 


< 5.6 x 10- 


31 


< 1.9 x 10 28 


< 6.2 x 10- 


-31 


< 1.7 x 10 28 


< 4.9 x 10~ 31 


< 1.4 x 10 28 




10 25 


< 1.7 x 10- 


23 


< 5.6 x 10 35 


< 3.4 x 10" 


-22 


< 9.1 x 10 36 


< 8.2 x 10" 23 


< 2.4 x 10 36 




1Q 20 


< 1.8 x 10" 


32 


< 5.9 x 10 26 


< 1.9 x 10" 


-32 


< 5.0 x 10 26 


2.9 x 10" 32 


8.5 x 10 26 


Extrapolated 


10 21 


< 1.8 x 10- 


32 


< 5.9 x 10 26 


< 1.9 x 10" 


-32 


< 5.0 x 10 26 


2.9 x 10~ 32 


8.5 x 10 26 


from 


10 22 


< 1.8 x 10- 


32 


< 6.0 x 10 26 


< 1.9 x 10" 


32 


< 5.1 x 10 26 


2.9 x 10" 32 


8.5 x 10 26 


2-10 keV 


10 23 


< 1.9 x 10- 


32 


< 6.4 x 10 26 


< 2.0 x 10- 


-32 


< 5.5 x 10 26 


3.1 x 10~ 32 


9.2 x 10 26 


Observed 


10 24 


< 3.5 x 10" 


32 


< 1.2 x 10 27 


< 3.9 x 10" 


-32 


< 1.1 x 10 27 


5.9 x lO" 32 


1.7 x 10 27 




10 25 


< 2.7 x 10" 


31 


< 8.9 x 10 27 


< 3.5 x lo- 


31 


< 9.5 x 10 27 


5.0 x 10" 31 


1.5 x 10 27 



Note. — Listed are the monochromatic fluxes and luminosities for the rest frame energy corresponding to 2 keV observed. The units 
are erg s — 1 cm -2 Hz^ 1 for flux, fx, and erg s — 1 Hz -1 for luminosity, lx ■ All fluxes and luminosities are upper limits, except for the 
3<r detection. 



RA=145.82335, DEC=54.28475, MJD = 52282, Plate= 769, Fiber=3+9 




z= 2.2309 V~, 0-000' (0.90),, Q50 

4000 saoa 6qod 7oao saoa saoo 

Wavelength [A] 



Fig. 1. — Observed optical spectrum of J0943+5417 from the Sloan Digital Sky Survey. 
The suppression of the spectrum here is a clear example of why it is difficult to correct the 
observed optical data. 



Table 4. Rest Frame 2500A Luminosity 



Object 


z 






/21.590A 


'~6900A 


^2500 from J 


^2500 from H 


(2500 from K s 








magnitude 


erg cm -2 s _1 Hz _1 


erg s — ^^Hz - 1 


erg s _1 Hz — 1 


erg s — 1 


erg s — ^^Hz -1 


J0943+5417 


2.23 


0.063 


14.251 


1.408X10 17 


1.6xl0 32 


9.1xl0 31 


1.3x10 s2 


1.6xl0 32 


J1352+4239 


2.04 


0.055 


13.861 


2.002X10 17 


2.0xl0 32 


1.5x10 s2 


1.5x10 s2 


1.8xl0 32 


J1723+5553 


2.11 


0.178 


14.097 


1.804X10 17 


1.9xl0 32 


7.8 x10 s1 


9.5 x10 s1 


1.6xl0 32 



Note. — Statistics for all three objects. The redshift, z, and galactic extinction, A u , are taken from Sloan Digital Sky Survey; K 3 is the 
observed magnitude from 2MASS; /21.590 is the observed flux in the K s band; '^gggoA * s * ne luminosity in the rest frame wavelengths of 
the K a band for each object; The last three columns contain the luminosity at 2,500A as extrapolated from the J, H, and K a bands. 
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333 normal quasars from ISteffen et al.l (120061 ). The sample was c hosen to cover as much 



FeLoBALs ( 


Brotherton et al. 


2005; 


Brotherton et al. 


( 


2005 


) assumed a 



galactic column density using ISchlegel et al.l ( 119981 ) and found lower than expected intrinsic 
values for their two FeLoBALs. They lie right on the bottom ed ge of and jus t below 



a 



(Aa m =0.3) the scatte r of expected a ox for the /2500 ~« ox relation from lSteffen et al.l (120 06). 



Rogerson et al.l ( 1201 ll ) compared the observed a ox of their FeLoBALs with the ISteffen et al. 
( 120061 ) study, and found the observed a ox to be ~ 0.48 below the expect ed scatter. From this 
they inferred an intrinsic column density. We follow the convention of ISteffen et al.l (120061 ) 
and use 2500A and 2 keV, although a m can be defi ned using other optical wavelengths 

In this section, we first examine 

calculated 



(12010( 1 . 



and X-ray energies, as shown by lYoung et al 

the dependency of a ox on the intrinsic absorption in our FeLoBALs, using a 
independently from the soft X-ray (0.2-2 keV, observed) and the hard X-ray (2-10 keV, 
observed) spectra. Then we analyze function of UV and X-ray luminosities, and 

draw conclusions about the intrinsic column density. 



4.1. Intrinsic Absorption and a ox 

We first study the relation between a ox and the intrinsic absorption in Figure |2J Each 
sub- figure corresponds to one FeLoBAL, and shows the entire input range for intrinsic column 
density. The values calculated from the 2-10 keV (hard) spectrum are depicted in black, and 
the 0.2-2 keV (soft) spectrum are depicted by red, filled symbols. Where monochromatic 
luminosities at 2 keV are upper limits due to the non-detection of a source, a ox is also only 
an upper limit, depicted in all Figures by arrows. Included are solid lines for the expected 



«ox~ ^2500 relation, as found by lSteffen et al.l ( 120061 ). and dashed lines to depict the expected 
scatter. 

At 2 keV the intrinsic luminosity and its uncertainty are dependent on the column 
density within the FeLoBAL itself, which causes intrinsic absorption that must be corrected 
for. We use a range 10 20 < Nh < 10 25 cm -2 for intrinsic column densities since the true 
quantity is unknown. This gives a different value of intrinsic luminosity, and therefore a ox , 
for each column density. The luminosities at 2500A are extrapolated from the rest frame 



optical luminosities, using a normal quasar SED from iRichards et al.l ( 120061 ). where we have 
corrected both the intrinsic extinction and Galactic extinction due to dust as discussed in 
§ [3j We find the uncertainty in the luminosities to be ~20%. 

In the two non-detected FeLoBALs (J0943+5417 and J1352+4239), for column densities 
up to 1 x 10 23 cm -2 , the soft limit is more stringent for a ox , although the difference between 
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the two limits diverges by < 0.2. At higher column densities, the hard spectrum is the 
appropriate limit to use. The divergence between the soft and hard limits is increasingly 
larger with increasing column density. J1723+5553 is a different story: only at Nh > 6 x 10 23 
cm -2 are the upper limits from the soft X-ray spectrum consistent with the 3a detection 
in the hard X-ray spectrum. The 3cr detection and its uncertainties for a ox calculated from 
the hard X-rays are >0.3 above the upper limits given by the soft X-ray energy range for 
column densities < 6xl0 23 cm -2 . Thus, the X-ray hardness ratio of J1723+5553 constrains 
the intrinsic column density to Nh > 6 x 10 23 cm~ 2 . 



4.2. Monochromatic Luminosities vs. a c 



Figure [3] shows the corr elation betwe e n l 2 flnp and a ox , with the best-fit linear regression 



for normal AGN found by ISteffen et al.l (120061) plotted. This relation is given by a ox = 
(-0.137 ± 0.008)log(/ 250 o) + (2.638 ± 0.240) Jsteffen et alJliooS ). The solid line shows the 
expected a ox , and the scatter is depicted by the dashed lines. Column density increases as 
a ox becomes less negative. The upper limits on a ox from the soft X-ray spectrum are shifted 
slightly to the right of the hard X-ray upper limits in Figure [3] for ease of viewing. The 
upper limits plotted for J0943+5417 and J1352+4239 are consistent with the normal AGN 



in the scatter of ±0.24 for N H > 4 x 10 6 and N H > 8 x 10 6 cm -2 , respectively. This range 
falls in the column density region (Nh > 1 x 10 23 ) where the hard spectrum provides the 
appropriate upper limit. For J1723+5553, combining the constraint from X-ray hardness 
ratio (Ah > 6 x 10 23 cm -2 ), the values of a ox fall within the expected scatter if the intrinsic 
absorption is in the range of 6 x 10 23 < Nh < 3 x 10 24 cm -2 . 

Figure H] shows the correlation between l^kev an d a ox , again with the best-fit linear 
regression a nd expected scatter, using the relation, a ox = (— 0.077±0.015)log(/2fc e v) + (0.492± 
0.387), from lSteffen et al.l (120061 ) . Column density increases from left to right, corresponding 
to rising intrinsic luminosity. For column densities Nh > 7 x 10 23 cm -2 for J0943+5417 
and Nh > 2 x 10 24 cm -2 for J1352+4239, a ox is consistent with the normal AGN sample. 
The upper limits given by the soft spectrum fall within (and above) the scatter in the 
high column density region, where the hard spectrum gives the better constraint. Our 3cr 
detection of J1723+5553 indicates that for column densities Nh < 8 x 10 24 cm -2 , J1723+5553 
is consistent with the SEDs of normal AGNs. 



Therefore, combining all constraints, J1723+5553 is consistent with normal AGN for 
column densities 6 x 10 23 < N H < 3 x 10 
are consistent with normal AGN for Nh > 7 x 10 
tively. 



:!4 cm" 2 , while a ox for J0943+5417 and J1352+4239 
23 cm -2 anc j j\j h > 2 x 10 24 cm -2 , respec- 
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Discussion 



5.1. Column Densities and Intrinsic X-ray Luminosities 

We find significantly high column densities for all three FeLoBALs studied, 6 x 10 23 < 
N H < 3 x 10 24 cm- 2 for J1723+5553, JV H > 7 x 10 23 cm" 2 for J0943+5417, and N H > 
2 x 10 24 cm -2 for J1352+4239. Only a small sample of FeLoBALs has been studied in X- 
rays previously. We summarize all the measurements in Table [5j which is divided into two 
subsamples: those that are radio loud, and those that are not. The radio loud objects tend 
to have lower column densities, which could be due to geometry or evolution. The X-ray 
emission may be linked to t he jets as wel l as fr om the corona right abo ve the central portion 



of the accretion disk (e.g., iMiller et al.l 120091 ) . The BAL wind (see iGallagher fc Everett 



(120071 ) for a detailed illustration), would then only absorb X-ray emission from the accretion 
disk region. We do not know for certain what complicated structure dependence for X- 
ray emission may exist in a radio loud FeLoBAL, but the column density results are still 
important and thus we include them for completeness. Our three FeLoBALs are radio quiet, 
and our constraints are consistent with the subsample that is not known to be radio loud. 
These objects in general are foun d to have A^h on the o rder of a few 10 23 cm -2 or higher. 



Notably, the most recent study by lRogerson et al.l (120111 ) found column densities for J2215- 
0045 and J0300+0048 to be N H ~ 10 24 , with J22 15-0045 just above the range we calculate 
for J1723+5553, but still consistent with our A^ h for the two non-det ections. The most recent 
calculation of column density for MrK 231 by iBraito et al.l (120041 ), using Bepp S AX, f ound 
A^h~ 2 x 10 24 cm -2 , in agreement with all three of our objects. I Green et al.l (120011 ) cite 



lower column densities (A^h > 6.5 x 10 22 cm 2 ), but they also noted that this is an average 
expected column density for a HiBAL, and that A% for a LoBAL would be higher. 

The fact that our 3a detection in J1723+5553 is in the observed 2-10 keV (rest frame 
6-30 keV) hard band rather than the soft X-ray band is significant, as this is expected from 
highly absorbed X-ray emission from the corona above the SMBH. The non-detection of 
the source in the observed 0.2-2 keV (rest frame 0.6-6 keV) band is consistent with our 
interpretation that we detect the direct X-ray emission in the observed 2-10 keV band, as 
high column densities will suppress the soft X-rays more than the hard X-rays. This is 
further supported by the fact that none of these objects is detected, even at the 3a level of 
significance, for the broadband observed 0.2-10 keV (rest frame 0.6-30 keV) energy range. 
If the observed 2-10 keV spectrum is reflected or scattered into our line of sight, the X-ray 
spectrum will be a flat power-law, and we would also detect the source in the observed 0.2-2 
keV band, which has a lar ger collecting area f o r Suzaku. Combining t he measurements from 
this paper and those from lBraito et al.l (120041 ); iRogerson et al.l (120111 ). it is most likely that 
the intrinsic Nr column density absorbing the X-ray emission in a radio-quiet FeLoBAL is 
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in the range from ~ 10 24 up to 10 25 cm 2 . This is about 1-2 orders of magnitude higher 
than the X-ray absorption range in BALQSOs (Njj ~ 10 22 — 10 24 cm -2 , e.g., Gallagher et 
al. 2006). 

Using a ox to provide limits on intrinsic column density also allows us to provide limits 
on the intrinsic X-ray luminosity. Using the combined constraints summarized at the end of 
Section |4T2| the intrinsic X-ray luminosity for J0943+5417 is /2koV — 9.8 x 10 26 erg s _1 Hz -1 ; 
for J1352+4239 it is / 2keV > 1.7 x 10 27 erg s" 1 Hz" 1 ; for J1723+5553 it is 4.4 x 10 27 < 



l okeM < 1.2 x 10 erg s Hz . Compared to the normal sample of AGN from lSteffen et al. 



( 120061 ). these luminosities fall into the middle to high range of X-ray luminosities for normal 
quasars. This is predicated on the fact that we assumed an intrinsic SED for a normal AGN. 
The broadband luminosities for these objects in the soft band are Lo.6-6 kev > 3.8 x 10 44 
erg s _1 , £ .6-6keV > 6.9 x 10 44 erg s _1 , and 1.8 x 10 45 < L .6-6keV < 4.9 x 10 45 erg s _1 , 
for J0943+5417, J1352+4239, and J1723+5553, respectively. The corresponding hard band 
luminosities are L 6 _ 3 okcV > 3.2 xlO 44 erg s" 1 , -^6-30 kev > 5.85 xlO 44 erg s -1 , and 1.5 xlO 45 < 
^e-so kev < 4.2 x 10 45 erg s" 1 . 



5.2. Thomson Scattering Implications 

Such high column densities are interesting physically, as they are either in or on the 
cusp of a regime where the Thomson scattering cross-section, <jt = 6.65 x 10~ 25 cm 2 , will be 
large enough to scatter incident photons. Although we probe the rest frame quasar spectra 
in the rest frame 6-30 keV (observed 2-10 keV) band, the Klein-Nishina correction to the 
Thomson scattering cross section is still negligible. For J1723+5553, using the expected a ox 
from the a oy rhkev relation, the optical depth is r ~ 1. For J0934+5417 and J1352+4239, 
the expected a ox value from the azox-hkeV relation gives an intrinsic Ah ~ 8.5 x 10 24 cm -2 , 
which indicates a Thomson scattering optical depth of r ~ 6. Such a large optical depth 
will also significantly modulate the UV flux (~400 times dimmer); ho wever, these o bjects 



are already more luminous in the UV than most of the AGN in the ISteffen et all (120061 ) 
study. As another check, we calculate the black hole mass of the quasars assuming they 
are emitting at 1/3 of the Eddington luminosity (e.g., Shankar et al. 2010). We extrapolate 
our rest frame luminosities to 5100A and assume a bolometric correction of 10.33 for 5100A 



( jRichards et all 120061 ) before using Equation to find the black hole masses. 



1 / L Edd \ 3L bo i 



l BH 



3 V 1.5 x 10 38 ^ 1.5 x 10 38 



M e * ; c S M (5) 



For the values reported here, all three FeLoBALs have Mbh ~ 1-4 x 1O 8 M . Even at the 
lower end of the column density constraints, Ah > 6 x 10 23 cm -2 , the optical depth would 
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Table 5. X-ray Observations of FeLoBALs 



FeLoBALs 



Object 


Detection 


Observatory 


N H (cm- 2 ) 


Reference 


J0943+5417 a 


no 


Suzaku 


> 7 x 10 23 


This paper 


J1352+4239 


no 


Suzaku 


> 2 x 10 24 


This paper 


J1723+5553 


no 


Suzaku 


6 x 10 23 - 3 x 10 24 


This paper 


SDSS J0300+0048 


no 


Chandra 


> 1.8 x 10 24 


Roeerson et al. (2011) 


SDSS J2215-0045 


no 


Chandra 


> 3.4 x 10 24 


Roscrson et al. (2011) 


Mrk 231 


yes 


BeppoSAX 


~ 2 x 10 24 


Braito et al. (2004) 


Q0059-2735 


no 


Chandra 


> 6.5 x 10 22b 


Green et al. (2001) 


Radio Loud FeLoBALs 


Object 


Detection 


Observatory 


N H (cm" 2 ) 


Reference 


SDSS J1556+3517 


yes 


Chandra 


< 9.6 x 10 23 


Kuncrt-Bairaszcwska et al. (2009) 


SDSS J2107-0620 


yes 


Chandra 


< 9.6 x 10 23 


Kunert-Bairaszcwska et al. (2009) 


SDSS J1044+3656 


yes 


Chandra 


< 9.6 x 10 23 


Kuncrt-Bairaszcwska et al. (2009) 


SDSS J0814+3647 


yes 


Chandra 




Miller et al. (2009) 


SDSS J2107-0620 


no 


XMM-Newton 


4 x 10 23 


Wane et al. (2008) 


FIRST 1044+3517 


yes 


Chandra 


~ 3 x 10 23 


Brotherton et al. (2005) 


FIRST 1556+4517 


yes 


Chandra 


~ 3 x 10 23 


Brotherton et al. (2005) 


FTM 1004+1229 


yes 


Chandra 


2.8 x 10 23 


Urrutia et al. (2005) 


FTM 1036+2828 c 


yes 


Chandra 


3.8 x 10 22 


Urrutia et al. (20051 


FTM 0830+3759 


yes 


Chandra 


2.7 x 10 22 


Urrutia et al. (20051 



Note. — Previous studies of FeLoBALs in the X-ray. We have added column density constraints when the study 
provided them. 

a Radio moderate. 

^ Green et al. I ibooil) note that this should represent an average HiBAL, and would be higher for a LoBAL. 
c Sometimes classified as a mini-BAL. 
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be t ~ .4, which would dim the intrinsic luminosity by one and a half times. At the column 
density corresponding to the expected a ox for J0943+5417 and J1352+4239, ~ 8.5 x 10 24 
cm -2 (from the otoyrhkev relation, and still within the limits for J 1723+5553), the optical 
depth of r ~ 6 would mean the intrinsic luminosity would be increased by 400 times, giving 
a black hole mass of Mbh ~ 3.2 x 1O 12 M . Such a large mass is not a physical possibility; 
therefore, it is unlikely that these objects are super-luminous in the UV, and we only see the 
tiny fraction of Thomson scattered emission. 

Since the X-ray emission region is expected to be smaller than the UV emission region 
from either variability (e.g., Chartas et al. 2001) or quasar microlensing (e.g., Dai et al. 
2003, 2010a; Pooley et al. 2007; Morgan et al. 2008; Chartas et al. 2009) arguments, it is 
instead possible that the X-ray absorbing material is located between the X-ray and UV 
emission region and only covers the X-ray e mitting region. This is consistent with the disk 
wind models of iGallagher fc Everettl ( 120071 ) and Murray et al. (1995), where an essential 
component of the model is the Compton thick shielding gas between the X-ray and UV 



emissi on, protecting the disk wind from being ionized by the X-ray emission. iRogerson et al. 



( 120111 ) also reached a similar conclusion, where the Thomson scattering optical depth was 
argued to have r > 3 (> 20 times dimmer). In this paper, the constraint is likely more 
stringent with r ~ 6 (~ 400 times dimmer). Aoki (2010) analyzed the UV spectrum of 
Jl 723+5553, and found the Ah column density in the UV wind was Ah > 5 x 10 17 cm -2 
using the curve of growth method with the unresolved Balmer absorption lines. Aoki (2010) 
pointed out that the covering fraction cannot be determined due to unresolved absorption 
lines, and assumed a covering fraction of 1. The column density would be higher than the 
given limit if the covering fraction is significantly less than 1. The limit derived by Aoki 
(2010) is seven orders of magnitude less than, but still consistent with, our X-ray column 
densities. If the X-ray and UV absorbing gasses are located in different regions, it can still 
be incorporated into both existing geometrical interpretations of FeLoBALs (e.g., Murray et 
al. 1995; Gallagher & Everett 2007) or evolutionary models (e.g., Fabian 1999). However, 
for evolutionary models the UV emission should come from the photosphere of the gas/dust 
cloud rather than the accretion disk. Such models have been simulated in recent studies 



(e.g., iCasebeer et al.ll2008l ). and support that FeLoBALs may be an evolutionary stage in 
the development of normal quasars. 

Another explanation is that FeLoBALs have intrinsic SEDs different from normal quasars. 
They could be extremely X-ra y weak compared to normal quasars, as in the case of the 
narrow- line quasar PHL 1811 ( iLeighly et al.l 120071 ). or the Narrow-Line Seyfert 1 galaxy 
WPVS 007 dGrupe et al]bo08h . This would nullify our A" H constraints for J0934+5417 and 
J1352+4239, since they are obtained by assuming a normal quasar SED. However, this would 
not explain J1723+5553, which we detect in the observed 2-10 keV band, but not in the ob- 
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served 0.2-2 keV band. If J1723+5553 were intrinsically X-ray weak, and X-ray unabsorbed, 
we would expect to see a flat-line power-law slope for quasars in the hard and the soft X-ray 
bands. This is not consistent with our observations. In addition, the Ah constraint from 
the X-ray hardness ratio is still valid for Ah > 6 x 10 23 cm -2 with the Thomson scattering 
optical depth r > 0.4, which could only marginally scatter the UV emission if the true Ah 
column density is close to the lower boundary. 

More X-ray detections of FeLoBALs are needed so that we are not restricted by using 
the broadband SED constraints. It is also possible that quasar host galaxies contribute 
a fraction of rest frame 4000-7000A emission; however, the contamination at this band is 
usually small. In addition, we correct the intrinsic dust extinction for FeLoBALs using the 
average dust extinction for LoBALs, which will underestimate the intrinsic flux and give less 
steep a ox values. A priori, it is unclear what the net effect of these two possibilities will be 
on the broadband flux; IR spectroscopy is needed to help resolve these issues. 



5.3. AGN Kinetic Feedback 

We calculate the kinetic feedback efficiency, = E^/LeoI, of the X-ray absorber in 
FeLoBALs using the location, luminosity, and Ah column density constraints obtained in 
this paper. In particular, is the kinetic feedback power, E k = Mv 2 /2 = 4ir[im p f c rN H v 3 
(e.g., Moe et al. 2009), where /z is the mean molecular weight, m p is the proton mass, f c , 
r, Ah and v are the covering fraction, location, column density, and velocity of the wind, 
respectively, and Lb i is the bolometric luminosity of the quasar. If the SEDs of FeLoBALs 
are consistent with those of normal quasars, we obtain a reasonable Ah ~ 3 x 10 24 cm -2 ; this 
is the upper end of the column density range for our 3a detection, and entirely consistent 
with our column density lower limits for the non-detections. We locate the wind between the 
UV and X-ray emission regions, ~ 40r 5 and r g = GM/c 2 , using the microlensing constraints 
of Dai et al. (2010a). For the covering fraction, we use the intrinsic FeLoBAL population 



range of 1.5% and 2.1%, depending on which model (iDai et al.ll2010bf ) is used, and assume 
the quasars are emitting at a typical 0.3 Eddington limit (Shankar et al. 2010). We are left 
with one uncertainty, the velocity of the X-ray absorber. The velocity of the BAL wind is 
usually measured in the UV wind, which can reach up to ~ 0.1c. Since the X-ray absorbing 
wind is located at smaller radii than the UV absorber, its velocity can be higher than the 
wind velocity measured from the UV spectrum. The only velocity measurement of the X-ray 
absorber is from the blue-shifted X-ray absorption lines detected in a few gravitationally 
lensed BALQSOs (Chartas et al. 2002, 2003, 2007), and this velocity can reach 0.3-0.8c. 
Since the X-ray absorption lines are mostly detected in mini-BALs, we could possibly be 
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looking through the edge of the wind, where the wind can be by fully accelerated. Thus, we 
consider the X-ray absorption line as providing an upper limit, and assume our wind velocity 
is 0.1-0. 3c, between the estimates from the two methods. We find the feedback efficiency, e^, 
for FeLoBALs is in the range (outflow velocity dependent) of either 0.2%-4.8% for a covering 
fraction of 1.5%, or 0.3%-6.9% for a covering fraction of 2.1%. To reach the minimum of 
5% required to explain the co-evolution between black holes and host galaxies (e.g., Silk & 
Rees 1998; Granato et al. 2004; Hopkins et al. 2005), with the other variables given here, 
we need a column density of Nn> 8 x 10 25 cm -2 or Nn> 6 x 10 25 cm -2 for low velocity 
wind with the different covering fractions. These are consistent with our lower limits on 
column density for the two non-detections, but out of reach for J1723+5553. However, for 
the high end of the conservative wind velocity range, only Nn> 3 x 10 24 and Nk> 2 x 10 24 are 
needed for f c =1.5% and 2.1%, respectively. Both of these column densities are consistent 
with the results for J1723+5553, and thus it is likely that FeLo BALs can contribute to 
AGN feedback. This becomes even more likely when models like iHopkins fc Elvisl (120101 ) 
are considered. They describe a "two-stage" feedback process, where a weak wind from the 
central engine energizes a hot, diffuse interstellar medium, which is then amplified when it 
it hits instabilities in a cold cloud within the host. This amplification outside of the central 
engine requires feedback efficiency as small as 0.5%. The BAL wind in FeLoBALs, therefore, 
is a promising candidate for the feedback process responsible for the co-evolution between 
black holes and host galaxies. More measurements of Nh are needed to further quantify 
how important the contributions are from FeLoBALs and BALQSOs to kinetic feedback 
efficiency. 
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Fig. 2. — Dependence of a ox on Intrinsic Absorption (Njj). Values calculated from the 2-10 
keV range are shown in black, while values calculated from the 0.2-2 keV range are shown 
in red filled symbols. The expected a ox from / 2 soo is depicted by the solid line, with the 
expected scatter depicted by the dashed lines. 
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